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We determine the behaviour of entanglement between regions of space in a Bose gas of fixed 
particle number around the critical temperature condensation. Long-range correlations develop 
in the Bose-Einstein condensate (BEC) phase transition and the aim here is to find out whether 
spatial coherence alone implies entanglement. We use a purity measure of entanglement to derive an 
entanglement witness that detects entanglement between two regions of space in the BEC. It is shown 
that spatial coherence between the two regions is necessary for entanglement with coherence and 
entanglement becoming equivalent only when the regions occupy the entire confining volume of the 
gas. The probabilities for Bosons to occupy the regions is the only other parameter that influences 
the amount of entanglement. We calculate explicitly the amount of entanglement between two 
regions for a cigar-shaped harmonic potential and find that it increases with decreasing temperature. 
A second entanglement witness is derived for entanglement between three regions of space, where 
again spatial coherence between each pair of regions is precursor to entanglement. It is shown that 
the state with the maximum amount of entanglement occurs when the probabilities for the Bosons 
to occupy the three regions are equal. 



I. INTRODUCTION 

The presence of entanglement in phase transitions has 
attracted considerable interest over the past few years. 
Understandably, quantum phase transitions were the 
natural place where one expected to find the non-local 
correlations of entanglement [TJ[H[3]. A quantum phase 
transition is the qualitative change in the ground state 
of a many-body system brought on at zero temperature 
by the change of an external parameter or coupling con- 
stant. At the quantum critical point long-range correla- 
tions develop in the ground state and an extra quantity 
known as the quantum order parameter is required to 
fully describe the system. The motivation for investigat- 
ing the entanglement in quantum phase transitions lies in 
the fact that the ground states of most condensed mat- 
ter many-bodied systems occupy huge Hilbert spaces and 
are difficult to construct. While one may not have full 
knowledge of such a ground state, one can still explore 
how the entanglement changes close to quantum critical 
points, thus accessing some aspects of the ground state's 
behaviour in critical regimes. Moreover, it has been ar- 
gued that entanglement can provide a physical picture 
and a deeper understanding of renormalisation methods 
f4^, quantum phase transitions [2] and their associated 
quantum order parameters [S]. 

Most research studying entanglement in phase transi- 
tions concerns discrete lattice type systems and to date 
hardly any research has focused on entanglement in the 
Bose-Einstein condensate (BEC) phase transition. Ini- 
tially, Simon [Hj demonstrated that any two regions of 
a uniform zero-temperature condensate of finite particle 
number are entangled. A recent paper [7 demonstrated 
that entanglement exists only below the critical tempera- 
ture for condensation, but a precise relationship between 
entanglement and the long-range order was not specified. 
The lack of research into entanglement in the BEC phase 
transition stems from the fact that there is no entangle- 
ment between individual particles in a non-interacting 



BEC. Usually entanglement is thought to exist between 
individual particles or their internal degrees of freedom, 
but in this case one has to instead consider entanglement 
between regions of space occupied by particles. 

The ground state of a non-interacting Bose gases is al- 
ready well understood, which means that we do not need 
to use entanglement to characterise it's behaviour. So 
why, in this case, is it useful to study entanglement in 
the BEC phase transition? While the BEC phase tran- 
sition is classical, i.e. it occurs by varying the temper- 
ature, quantum effects develop on a macroscopic scale 
without the presence of interactions. This makes a non- 
interacting BEC a remarkably simple platform for estab- 
lishing concrete connections between entanglement and 
the long-range correlations that develop in phase tran- 
sitions. If we can fully understand how entanglement 
emerges in the simple setting of the BEC phase transi- 
tion, we can use these results as a guide for understanding 
how entanglement and long-range correlations are related 
in more complicated systems. For instance, Osborne and 
Nielsen [U|3] conjecture that entanglement is responsible 
for the long-range correlations that develop in the trans- 
verse Ising and XY models at the quantum critical point 
and that at such critical points the system is maximally, 
or critically, entangled. 

In this paper we will develop a thorough understand- 
ing of entanglement in the BEC phase transition and as 
a consequence we will be able to provide evidence for and 
against such statements. We will see that spatial coher- 
ence of a Bose gas is necessary for entanglement between 
regions of space. The probabilities for the Bosons to oc- 
cupy the regions is the other crucial parameter affecting 
the entanglement. So while there may be spatial coher- 
ence between two points in space, these points will not 
be entangled because there is almost a zero probability 
of Bosons populating the points. In this way we will be 
able to understand what portion of the long-range corre- 
lations are actually entanglement. 

The remainder of this paper is organised as follows. 
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We will first discuss the concept of long-range order and 



spatial coherence of a BEC in section I A followed by 
brief discussions of spatial entanglement in section |IB| 
and an overview of an entanglement measure based on 
the purities of states in section |I C[ In section [ll] we will 
calculate the entanglement between two regions of space 
around the critical temperature for condensation. An in- 
equality relating the entanglement to spatial coherence 



will be derived in section III We will compare our the- 



oretical results for entanglement to recent experimental 
results measuring the spatial coherence of a BEC in sec- 
tion |IV[ To fully characterise entanglement in the BEC 
phase transition we will consider multi-mode entangle- 
ment in section |V| Finally in |VI| we will conclude by dis- 
cussing the results of this paper in the context of previous 
results concerning the role of entanglement in quantum 
phase transitions of discrete systems. 



A. The BEC phase transition 

A simple but important example of a phase transition 
is the formation of a BEC. As a Bose gas is cooled close 
to absolute zero, it passes a critical temperature where 
the chemical potential becomes comparable to the ground 
state energy causing a macroscopic number of Bosons to 
collapse into the ground state of the confining potential. 
As a result, all of the Bosons behave as a single wave 
and the phase of the condensate becomes locked. The 
wavefunction of the condensate then plays the role of the 
order parameter. The appearance of an order parame- 
ter signifies that long-range correlations have developed 
between two distant points of a condensate. 

Such long-range correlations were described mathe- 
matically by Penrose and Onsager [S] in terms of the 
long-range order (LRO). LRO is present in a condensate 
when the one-body reduced density matrix 



pi(f,r") = tr[^'T(r)^'(r')p] 



(1) 



is finite as |r— r '| oo, where ^''^(f^ and ^'(f^ are the cre- 
ation and annihilation operators of particles at point r. 
This definition holds for interacting and non-interacting 
gases, but because pi{f,f') tends to the ground state 
density, p\{f,f') no, for |r — r'| oo, its value is 
diminished by particle interactions. The one-bodied re- 
duced density matrix, r'), measures the presence of 
off-diagonal matrix elements in the state and determines 
the quantum coherences between two points. The aim of 
this paper is to investigate the how the spatial coherence, 
described by pi(r, r'), are related to entanglement. 

Let us briefly consider whether the one-bodied re- 
duced density matrix, p(r,r'), is a good measure of 
entanglement. If we imagine one particle distributed 
at two points, r and r ', with equal probability, the 
state of the two points occupied by particles would be 
IV'i) = 1/V2(|10) + |01)), where |01) = |0)a ® \1)b- A 
two particle state is IV'2) = l/%/3(|20) + |11) + |02)) and 
a N particle state would be IV'jv) = 1/\/n{\NG) -I- - 



1, 1) -I- . . . -I- |1, iV - 1) -I- |0A^)). The N particle state, 
\iPn) , clearly contains more entanglement than the one or 
two particle state, but the normalised one-bodied reduced 
density matrix is identical for all three cases. So while 
the one-body density matrix can detect entanglement, it 
cannot differentiate between entanglement formed by one 
or N particles occupying the regions. 



B. Entanglement between spatial modes 

What does it mean for two regions of space to be entan- 
gled? We know from the simple form of Hamiltonian of a 
non-interacting BEC, H = ^^E^a^Ji-^, that individual 
particles are decoupled in momentum representation and 
are therefore not entangled to one another. However, one 
could also view this Hamiltonian in position representa- 



tion. Here the operators, al, a^r, that create and destroy 
particles in mode fc, transform as at — J 

where <j>j^{x) is the eigenvector for mode k and tp^{x) cre- 
ates a particle at position x. In position representation, 
one can split the integral into multiple regions so that 
there are pairwise couplings between all regions of space, 
H = E.,, /, dr df cl^f:{r)<j>l{ni'Hmr')- This 

is in fact what leads to entanglement between regions of 
space. 

Spatial entanglement therefore takes the form of non- 
local particle number correlations between regions of 
space or spatial modes. The spatial modes are entan- 
gled and not the particles themselves. In order to ensure 
the correct Hilbert space structure one must work with 
second quantised modes and not in the first quantised 
picture. In second quantisation the Hilbert space of three 
spatial modes occupied by Bosons have a tensor product 
structure. One is then able to trace out each mode to 
obtain a state of just two regions of the gas. For a more 
detailed discussion of this and other aspects of spatial 
entanglement see references [3 [12 [T31 [TH [TS] . 



C. Entanglement of continuous variable states 
using purity measurements 

In this paper we will use a measure of entanglement 
based upon the purities of the global and reduced states 
of the spatial modes. The purity, /i, measures the degree 
of information contained in a quantum state, p, and is 
defined as p = tr[p^]. The purity is related to the mixed- 
ness, Sl{p), of the state by Sl{p) = [N/{N - p), 

where N is the number of modes. The purity ranges 
from one for pure states to 1/D for fully mixed states, 
1/D < p < 1, where D = dimTi, is the dimension of 
the Hilbert space. For a state, pab, of two modes A and 
B, the global purity is defined a,s p — tr[p^g] and the 
reduced purity of the ith mode is pi = tr[pf], where pi is 
the reduced density operator of mode i = A or B. 
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It is well known that the concepts of entanglement 
and the information encoded in a quantum states are 
closely entwined. For instance, in pure states, entangle- 
ment and the lack of information, or mixedness, of the 
reduced state of a subsystem are equivalent. A maxi- 
mally entangled pure state will have maximally mixed 
reduced states. Thermal states on the other hand have 
a mixed global state. Entanglement therefore exists in 
mixed states if the global mixedness of a quantum state 
is less than some function of the mixedness of the reduced 
states of the subsystems. 

An entanglement measure in terms of global and re- 
duced purities for continuous variable Gaussian states 
[HI is 

When this inequality is satisfied, there is entanglement 
between modes A and B. The upper bound represents 
a physicality condition on the quantum state, pad- Just 
below the lower bound there is a small region where en- 
tangled and separable states coexist. We shall see shortly 
that for systems with a large number of non-interacting 
particles, the region where separable and entangled states 
both exist goes to zero. 

In what follows we shall choose the spatial modes so 
that the reduced purities of A and B are equal pA — 
Pb — Pr- The entanglement criterion then reduces to 



at as 11,-) = a- 10). But, as entanglement does not 

k ^ k 

exist in momentum representation, we transform into 
position representation as jl^:) = J dr(j)j:{'r)tp^{'r)\0), 
where the operator V'^CO creates a particle at posi- 
tion f. To form the three spatial modes A, B and 
C, the integral in jl^:) is split into the three parts 

giving |lfe)^= {J^dr(f>j:{T^i!^r) + Jgdr(f>j:{r)ti;^r) + 
J^dr(j)j:{r)^lj^r))\0). The state of a single particle is 
p = -ffel-ifc)(-ifcli where each momentum mode is 
weighted by the probability of finding it in that mode, 
Pj: = n^/N, where J^k^k ~ ^ ^^"^ '*fc ^s the Bose distri- 
bution [IB]. 

So far we have found the state of one particle expressed 
in terms of three spatial modes. To arrive at pab, re- 
gion C is traced over. We can perform the trace on each 
momentum mode separately, which results in the state 

— trc[p^*^^] for mode k, where = |l^:)(lj:|. Ex- 

(k) 

plicitly p\^ is 

PaI= E [dr f dr^M^rkin^^r^mmi^ir^) 

i,]=A,B-'^ •'3 



Pr 



<P<1, 



(3) 



and will be used shortly to derive our entanglement wit- 
ness. 



II. ENTANGLEMENT BETWEEN TWO 
SPATIAL MODES AROUND THE CRITICAL 
TEMPERATURE FOR CONDENSATION 

We will now calculate the entanglement between two 
spatial modes in a BEG. Let us consider a non-interacting 
Bose gas of fixed particle number, N , at temperature T, 
placed in some confining volume V{r), where the parti- 
cles occupy energy levels determined by the eigenfunction 
4>kif)- The results that we present here are for any con- 
fining potential; we will present a specific example in sec- 
IV The purity measure, d3]), requires knowledge of 



tion 

the states of one and two spatial modes, pA and pab re- 
spectively, where regions A and B do not together occupy 
the entire volume V{r). To arrive at, pA and pab, the 
state of the gas is mathematically divided into three re- 
gions A, B and C, as this allows us to speak about spatial 
modes. In reality the gas remains unchanged as we are 
not considering physical partitions at all. Region C will 
then be traced out leaving the state, pab = ^t:c[pabc]i 
of regions A and B that may be entangled. 

To find the state, pab, of one particle, we first cre- 
ate a particle in mode k using the creation operator 



df|0fc(r)noO)(OO|, 



(4) 



c 



where the vector |00) = |0)a ^ |0)b is the vacuum for re- 
gions A and B. The final state, pab, is therefore pab — 
J2k -^kPAB ^^'^ likewise for pA- Let us take a closer look 
at how one performs the trace in this situation. The op- 
erators act on point r which is located within 
one of the three spatial modes. The trace of the ith spa- 
tial mode acts only on operators whose position vector r 
is located in mode i. One can write the trace of the ith 
mode as tr,[...] = E„=o,i /, rfr(0|(^(r))" ... (^t(^-))n|o). 

As we are considering a gas of non-interacting Bosons 
with Hamiltonian H = "^k -^k^tcij:, it is straightfor- 
ward to generalise to the case of N particles. We can 
write the N particle state in terms of the single particle 
state as p'-^^ — p^-^ . To prove that this is equivalent 
to the state for BEGs found in text books consider one 
Boson in first quantisation. One can place N particles 
into this state, ensure it is correctly symmetrised and 
then move into second quantisation. Equally one could 
place N particles in the first quantised state and move 
directly into second quantisation, which results in the 
state above. For TV particles, the entanglement criterion 
^ can be expressed differently in terms of the purities 
for the single particle density matrices, p = tr[pAS^] and 
Pr = tr[pA^]- The purity, pN of N non-interacting par- 
ticles is pn = tr[{p^'^)'^] = p^ , where p is the single 
particle purity. Using this relationship the inequality for 
entanglement becomes p — Pr > for large N. If the 
global purity is larger than the reduced purity, the two 
spatial modes A and B are entangled. Indeed the upper 
bound, /i < 1, is always satisfied if the state, p, is physi- 
cal to begin with and for large N {N > 1000), separable 
and entangled states never coexist. 
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The global and reduced purities for a single particle 
can now be evaluated. The global purity, fi = tr[pAB'^], 
is 



III. THE RELATIONSHIP BETWEEN 
ENTANGLEMENT AND SPATIAL COHERENCE 
IN A EEC 



i=A,B,C 



df 



df' p^{f, f), 



(5) 



where pi = J,^dr\(j)j:{f)\'^ is the probability 

for occupying the ith spatial mode and p(f,f) — 
Sfc -ffc '?^fc(^<^^('^)- It is immediately clear that the state 
PAB is not pure, as tracing out region C destroys the co- 
herences between it and regions A and B. Equally, the 
reduced purity for the state of one region contains no 
coherence terms and for region A reads 



PA = tr[p\] = {pB +Pcf +Pa = p(p - 1) + 2 (6) 



Note that for simplicity the reduced purities for A and 
B are chosen to be equal, pa = l^B — Mr- 

Using the entanglement criterion, p — p^ > 0, we can 
write down a quantity, £, that will be positive when there 
is entanglement between spatial modes A and B and is 
negative when the two modes are separable. The quan- 
tity £ can therefore be thought of as an entanglement 
witness and is given by 



£ ^ -pa - 2p) 



df 



df p'^{f,f). 



(7) 



When £ > 0, the two regions A and B are entangled. 
Eq. is a general expression for non-interacting Bose 
gases as one can consider a gas with any normalised mo- 
mentum distribution Pg, trapped in any geometry in any 
dimension. The second term, or coherence term, in £ 
measures the spatial coherence between regions A and 
B. The first term in £ is parameterised by the probabil- 
ity, p, to occupy region ^ or i? and takes into account 
the fact that while there may be coherences between two 
regions of any size, there is not necessarily entanglement. 
For entanglement to exist the purity of the global state 
has to be greater than the purity of the reduced state. 

If the regions A and B are half of the volume of the gas 
each, the probability to be in each region is p = 1/2 and 
the second term in £ is zero. In this case entanglement is 
equivalent to spatial coherence between the two regions. 
This is particularly interesting as a small amount of co- 
herence, and hence entanglement, will remain between 
two such regions above the condensation temperature. 
On the other hand, if we choose two very small regions 
while there will be spatial coherence between them, there 
will be a relatively low probability for the Bosons to oc- 
cupy them and there will be no entanglement. We will 



The normalised one-body reduced density matrix, 
Pi{r,r') = Eki^k/N) (t)j:{f)(j)t{f) is identical to 
p(f,f) — 'E%Pi:4>%{'r)4'*j^{f') used in the expression for 
entanglement £, Eq. ([t]). The entanglement witness, £, 
is already expressed in the language of LRO and spatial 
coherence. It is evident that spatial coherence is neces- 
sary for entanglement between regions A and B, as the 
second term in £ is, of course, always positive. 

Rather than looking at the coherences between two 
points, it may be more appropriate to define a quantity, 
Ou = tr[V']i'/'B p] that measures the coherence between 
the regions A and B. The operators 7/;^(r) and '4>{f) are 
averaged over a detector profile g(f)in the regions, result- 
ing in ~ J^dfg{f)'4!^{f) and ipB = /g '^'^5*(^V'(^ 
that create and destroy particles in regions A or B. The 
quantity Od is given by 



Or 



tr[-ip\ifiBp] 



df / df'g{f)g*{f')p{f,f'). (8) 



From elementary statistics we know from the positiv- 
ity of the variance, a'^ — {X"^) — (X)^, that the expec- 
tation value of some quantity squared is greater than 
the square of the expectation value of the quantity, X, 
namely {X'^) > {X)"^ . Thus, the quantity Od squared 
is always less than the first term in £- the expectation 
value of the one-body reduced density matrix squared. 



df / df'p\{f,f') > O 



(9) 



The left hand side of this inequality is the coherence term 
in the entanglement witness, £, which is always greater 
than or equal to the newly defined coherence term Od- 
Using this fact we can define a new quantity £o,j that 
also acts as an entanglement witness. 



£on=0l-p{\-2p). 



(10) 



plot £ for a specific system in section IV 



If the new quantity, £odi i^ greater than zero, £0^ > 0, 
the previous entanglement inequality, £ > 0, is automat- 
ically satisfied and there is entanglement between regions 
A and B. It is worth expressing the entanglement witness 
in terms of Od as this quantity only depends upon the 
one-body density matrix, p{f, f), rather than its square, 
and is straightforward to measure experimentally, see for 
instance [8]. 



IV. A SPECIFIC EXAMPLE: EEC IN AN 
HARMONIC TRAP AND COMPARISONS WITH 
EXPERIMENT 

To bring clarity to the results presented in the previ- 
ous sections we now present a specific example. We will 
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FIG. 1: The amount of entanglement as measured by the av- 
erage logarithmic negativity [17] between two spatial modes 
in a harmonically trapped non-interacting BEG of fixed par- 
ticle number above and below the transition temperature for 
condensation [Tc = 440 nK) . See text for the trapping param- 
eters. The amount of entanglement decreases with increasing 
temperature of the gas. Due to the Gaussian ground state 
of the harmonic oscillator the amount of entanglement also 
decreases as the two regions move further apart. 



consider a typical experimental set up, namely a BEC 
confined in a harmonic cigar shaped trap. The exact 
trapping parameters that we use here have been taken 
from an experimental paper that measures the spatial co- 
herence of a BEC |8] as this will allow us to make direct 
comparisons with experiment. In this case the transition 
temperature for condensation is Tc « 440nK, the trap- 
ping frequencies are uoax = 27rx 13 Hz and tOper = 27rx 140 
and the number of Bosons in the gas is « 4 x 10^. We 
take a one-dimensional slice through the trap, define re- 
gions A and B of size La — Lb ~ 350nm and trace out 
the rest of the gas. The amount of entanglement as deter- 
mined by d?! and quantified by the average logarithmic 
negativity [17^ can be seen in Fig. [l] Entanglement does 
not exist above the critical temperature for condensation 
and the amount of entanglement between the two regions 
increases as the temperature decreases. There is no en- 
tanglement between the regions when they are separated 
by more than Az — 600nm. 

However the harmonic nature of the trap hides some of 
the interesting features of the entanglement. Therefore 
let us review the results presented in Simon A BEC 
was considered at zero temperature, where the Bosons 
occupied all points with equal probability. The same 
amount of entanglement was found to exist between two 
regions of a given size at any separation, i.e. two re- 
gions that were touching or two regions that were very 
far apart have the same amount of entanglement between 
them. This is because in position representation there are 
pairwise couplings between all regions of space (see sec- 
tion I B ) and in [6^ the probability to occupy any region 



> 



1.0 ■ 
0.8 ■ 
0.6 - 
0.4 
0.2 
0.0 



* + + 



• Az = 185 nm 

■» Az = 325 nm 
□ Az = 600 nm 



» * 4 







100 200 



300 400 
r(nK) 



500 600 700 



FIG. 2: Spatial coherence of a BEC above and below the 
transition temperature Tc as measured by Bloch in et al. 
[S]. The temperature dependency of the measured visibility 
is plotted for three separations. Above Tc (=430nK) the visi- 
bility is zero for slit separations 325nm and 600nm. A sudden 
increase in the visibility occurs below Tc with the maximum 
visibility being attained close to T = 0. The visibility does 
not vanish above Tc for the slit width Az = 185nm as the 
coherence length at the transition temperature is still larger 
than the slit width. 



of space was equal. Here entanglement decreases with 
increasing region separation due to the Gaussian ground 
state wavefunction. As the regions move further apart 
the probability for Bosons to occupy the regions drops, 
and the state of regions A and B becomes increasingly 
mixed. 

We can compare our results for spatial coherence to 
the experimental results. Let us start by describing in 
more detail the experiment by Hansch et al. ^ that 
measures the spatial coherence, pi(r, r'), of a weakly in- 
teracting Bose gas. In this double slit type experiment 
two matter waves were extracted from two regions in the 
trap separated by Az. The matter waves were allowed 
to interfere and the degree of spatial cohrence was deter- 
mined by the the visibility of the fringes. The position 
of each 'slit' region was determined by the frequency of 
r.f. field that created transitions between magnetically 
trapped and untrapped states. The width of the slits 
were small, as the r.f. field only coupled weakly to the 
gas. The results are given in Fig. [2] where the spatial 
coherence is plotted for separations up to Az — 600nm. 

We can make comparisons between the experimental 
results measuring spatial coherence between two regions 
m Fig. [2] and our theoretical results for entanglement 
and spatial coherence in Fig. [3j Let us first take a look 
at the spatial coherence in both cases. Spatial coher- 
ence defined hy Od Eq- [S) is denoted by the solid lines 
in Fig. |3] There is greater spatial coherence between 
two regions that are closer together. In our theoretical 
results the spatial coherence for all three separations is 
very similar. If we compare this to the experimental re- 
sults of Fig. [2] we see that the spatial coherence for sep- 
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FIG. 3: The entanglement and spatial coherence at three sep- 
arations in a non-interacting BEG. The amount of entangle- 
ment between two spatial modes of size La = Lb = 350nm 
is denoted by the dotted lines for slit separations of Az = 
185nm, Az = 325nm and Az = 600nm. The spatial coherence 
(LRO on legend) as defined by Od , Eq. ([sjl also for the three 
same slit separations is is denoted by the solid lines. Spatial 
coherence is always present when entanglement exists. The 
amount of entanglement depends strongly on the separation, 
and hence population, of the spatial modes as the spatial co- 
herence does not change much between the Az — 185nm and 
Az = 600nm. There is no entanglement for the slit separa- 
tion Az = 600nm because, while there is spatial coherence, 
the global purity of the two regions is lower than the reduced 
purity, due to the increased separation of the regions. Golour 
online. 



arations Az = 600nm and Az — 325nni decreases much 
faster than our predictions. This is because Bosons in a 
real gas interact with each other, which depletes ground 
state energy level and hence reduces the spatial coher- 
ence. We have not included interactions in our analy- 
sis as they only reduce the amount of spatial coherence 
by a small fraction, which would in turn only effect the 
amount of entanglement present by a small amount. 

The entanglement in Fig. [3] is represented by dashed 
lines and exists only when spatial coherence exists. Ini- 
tially the amount of entanglement between the regions 
is governed by the behaviour of the one-bodied reduced 
density matrix, p{r,r). When the spatial coherence, 
p(r, r), between two regions falls below the probability for 
the Bosons to not occupy the regions the entanglement 
falls sharply to zero. The global purity is larger than the 
reduced purity to begin with due to the presence of the 
coherence terms, but as the gas is heated the coherence 
terms reduce (solid lines Fig. [3| until the global purity 
and reduced purity are equal. The entanglement then 
disappears. There is no entanglement between the re- 
gions when their separation is greater than Az = 600nm, 
even when spatial coherence between the regions is max- 



imum (at T = 0). This is due to the Bosons gathering 
around the centre of the trap so that the probability to 
occupy the regions is very low. Here the state of two re- 
gions is very mixed, even at T = 0. On the other hand if 
the regions A and B occupy the entire trap, the amount 
of entanglement between them depends only on the spa- 
tial coherence between A and B and decreases steadily 
to zero as the spatial coherence goes to zero. 



V. MULTI-MODE ENTANGLEMENT IN A BEG 

One must carefully choose the set of spatial modes to 
ensure they detect the entanglement in the gas. We argue 
now why a BEG at zero temperature is entangled, even 
though our choice of spatial modes may not detect the 
entanglement, such as two regions separated by Az = 
600nm. 

At zero-temperature, T — 0, all of the Bosons occupy 
the ground state and the total state of the system, |^), 
is a pure state. If the state, |^), were separable with re- 
spect to three spatial modes it would be a product state 
as follows l^*) = \tpA) ® IV'b) ® IV'c)- If region C were 
then traced out, the resulting state, I^Pab), would also 
be a pure state, I^'ab) = IV'a) <8) lips) and would conse- 
quently have unit purity, /x = tr[yO^^] — 1. However in a 
BEG, the purity, fi, of the state, pab, is always less than 
one, p — tr[p\g] < 1, unless regions A and B together 
occupy the entire confining volume. Therefore, when the 
purity is less than one, p = tr[p^g] < 1, we can con- 
clude that the state of three spatial modes, j^*), was not 
a product state to begin with. This means that a BEG of 
fixed particle number at zero-temperature is always en- 
tangled, because a state of separable regions would have 
unit purity. We can also investigate multi-mode spatial 
entanglement to fully understand how the choice of spa- 
tial modes affects the detection of entanglement in the 
BEG phase transition. 

We will investigate entanglement between three spa- 
tial modes. A, B and C, that together occupy an en- 
tire confining volume, V, at finite temperatures, T. We 
can no longer use the purity measure of entanglement as 
this is only valid for bipartite states. However, there are 
suitable entanglement witnesses that detect genuine tri- 
partite entanglement. An entanglement witness, Hw is 
an observable that, by our notation, has a positive ex- 
pectation value for entangled states, trpv^PASc] > 0; 
and a negative expectation value for separable states. 
We will use the witness Uw — \W3){W3\ — 2/3, where 
iWs) = 1/V3(|100) + |010) + |001)), as our single particle 
state, PABC, belongs to the W class of tripartite entan- 
gled states. As we are studying a continuous variable 
system, we shall write the entanglement witness, Hwi in 
terms of field operators, -(/'^(r) and •ip{r), as 



E df dr'gir)g*{r')4,\r)\0){0\4,{f)-2, 



(11) 
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FIG. 4: Tripartite entanglement, as detected by the entan- 
glement witness W, of a Bosonic gas trapped at zero tem- 
perature, r = 0, in an infinite square well, as a function of 
the lengths. La, Lb and Lc, of the spatial modes A, B and 
C respectively. The lengths are given as fractions of the well 
length, L, and the legend indicates the length of the spatial 
mode A, La- Entanglement exists when the probability to 
occupy the regions are of the same order. Entanglement does 
not exist between the three spatial modes. A, B and C, when 
any one region is larger than 2/31/. Region B can become 
much smaller. Lb ~ 0.02 L, than the other two regions and 
entanglement still exists. See text for discussion. Colour on- 
line. 



where we have used = ^a\0) = J^drg{r)il)^r)\0). 
Here g{r) takes the role of a detector profile that specifies 
how we average over the set of points included in the 
mode. 

We are now in a position to construct an entanglement 
witness, W = trpv^PABcli that detects tripartite entan- 
glement of a single-particle in the BEC between the three 
modes. The witness, W, is 



^ I dr I dr' g{r)g{r')pi{r,r') 



(12) 



ij=A,B,C 



where pi{r,r') = Y^k^kh'^^'t^li^')- When W > 
there is tripartite entanglement between the three spa- 
tial modes A, B and C. As we would like to know how 
entanglement is related to spatial coherence and the one- 
bodied reduced density matrix pi(r, r '), it suffices to cal- 
culate the entanglement of a single particle in the Bose 
gas. This is because while pi{r, r') detects entanglement, 
it does not differentiate between single-particle and A^- 
particle entangled states, (see section lAl. 

The entanglement between three spatial modes. A, B 
and C, depends solely on the behaviour of the spatial 
coherence between each pair of regions, as here we have 
not traced out any of the gas. If we consider a state at 
very high temperature, T >> 1 every point in space will 
be separable from the points immediately next to it, as 
the de Broglie wavelengths of the atoms become infinites- 
imally small. Correspondingly, all of the off-diagonal co- 



FIG. 5: The maximum temperature for entanglement, as de- 
tected by the entanglement witness W, between three spatial 
modes of a Bosonic gas in a three-dimensional uniform poten- 
tial of volume y = L^, as a function of the size of region B, 
Lb ■ Note that region B is situated between regions A and C, 
whose lengths are taken to be equal, La = Lc- The highest 
temperature for entanglement is Tmax = 196nK, which is be- 
low the critical temperature for condensation, Tc ~ 270nK, 
which occurs when the probability to occupy all three regions 
is equal. See text for further discussion. 



herence terms present in the witness W disappear. The 
three regions are of course separable, as W — —1. 



A. Zero temperature tripartite entanglement 

On the other hand, at zero temperature, all of the 
Bosons occupy the ground state and coherences stretch 
across the entire system. Whether the three spatial 
modes are entangled in this case depends only on the rel- 
ative sizes of the three regions L^, Lb and Lc- We inves- 
tigate the entanglement between three spatial modes as 
detected by our entanglement witness, W, as a function 
of the region sizes. We consider again a Bosonic gas of 
fixed particle number, N, trapped at zero-temperature, 
T = 0, in the ground state, (f>o{r) — yj2 / L sm{T:x / L) , of 
a one dimensional infinite square well of length, L- 

Entanglement exists only when the probabilities to oc- 
cupy the regions are of the same order. See Fig. |4] Re- 
gions A and C are located on the very left and the very 
right of the well respectively and region B sits between 
them. Entanglement exists when the sizes of regions A 
and C are between 0.65 L < La,c ^ 0.125 L. Region 
B, however, can become much smaller. Lb ~ 0.02 L and 
still be entangled to both regions A and C, due to it's 
position in between regions A and C. This is because the 
probability for particles to occupy a small region in the 
centre of the trap is much higher than the probability to 
occupy the same small region at the edge of the trap. 

To see more clearly why this is the case let us consider 
a general W state, \^l^w) = a|100) /3|010) 7|001). Let 
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the coefficient, a become much less than the other two co- 
efficients, a << /3,7. The first term, whose coefficient is 
a, can be set to zero. The first subsystem (denoted by the 
first position in the kets) is then no longer entanged to the 
other two and the M^- witness, Ilw, no longer detects any 
tripartite entanglement. This is what, in fact, happens 
here when one or two of the regions becomes relatively 
small compared to the remaining one(s). However, the 
ground state of the gas is still entangled even when our 
witness W does not detect it. This is because the spatial 
modes chosen were not the optimum choice, i.e. we could 
have uncovered entanglement between smaller regions of 
space by choosing multiple (n > 3) small spatial modes. 
It is an open problem how the amount of entanglement 
scales with an increasing number of regions. 



B. Thermal tripartite entanglement 

In between high and zero temperatures there is an in- 
crease in the thermal de-Broglie wavelength of the atoms. 
The atoms spread over neighbouring regions of space and 
quantum coherences form between them. At a certain 
temperature the coherences between the three spatial 
modes will become big enough so that the regions be- 
come entangled. We will see that entanglement forms 
below the critical temperature for condensation. Let us 
consider a non-interacting Bose gas of = 30, 000 atoms, 
confined in a three-dimensional uniform potential of vol- 
ume, V = L^, where L = SlSnm. The space is divided 
into three regions, A, C and B of lengths La = Lc and 
Lb respectively along one axis and of length L along 
the other two axes. Region B is situated in between re- 
gions A and C. The atoms occupy the momentum modes 
(t>-j:{'r) = yj 2 / L sixiikr) . The critical temperature for con- 
densation is Tc — 270nK. 

Fig. [5] shows the maximum temperature, Tmax, for 
which entanglement exists as a function of the size of 
region B, Lb (which is also the separation of regions 
A and C), as determined by the entanglement witness, 
W, Eq. ( [T2| ). The warmest temperature that entangle- 
ment survives at is Tmax = 196nK, which corresponds to 
the region sizes La — Lc — 0.4 L and Lb — 0.2 L. At 
this point there is approximately equal probability, P, 
for the Bosons to occupy any region, i.e. Pa = Pc ~ Pb- 
At zero temperature this spatial configuration of modes. 
Pa ~ Pc ~ Pb, represents a state with the maximum 
amount of possible entanglement, as maximally entan- 
gled states usually posses equal probability amplitudes. 
I.e. the maximally entangled Bell states all have prob- 
ability amplitudes a — 1 j \pl. As the temperature in- 
creases from zero Kelvin, a fraction of the Bosons occupy 
higher momentum modes and the entangled state of three 
spatial modes becomes increasingly mixed. At some tem- 
perature, these quantum states become so mixed that 
they are no longer entangled. States with only a small 
amount of entanglement at T = need only a small 
amount of extra mixedness to become separable so are de- 



stroyed by only a small amount of heat. In contrast con- 
sider the entangled states with some maximum amount 
of entanglement at T — 0. The entanglement in these 
states survives to higher temperatures as they tolerate 
more mixing before they become separable. 

As the number of Bosons, N , in the gas increases the 
critical temperature for condensation, Tc also increases 
for a fixed volume. How does the maximum temperature 
for single-particle entanglement, Tmax, change as N in- 
creases? Naturally, Tmax increases with N and Tc, but 
the ratio, T^ax/Tc also tends to unity, Tmax/Tc 1. 
Thus, as the number of Bosons in the condensate, N , 
becomes very large, N — > 00, the maximum temperature 
for entanglement Tmax tends to the critical temperature 
for condensation, T,-aax — > Tc- In this limit, and for the 
optimum region configuration described above i.e. when 
Pa = Pb — Pc, spatial entanglement and spatial coher- 
ence become equivalent. In other words entanglement 
between three regions of space and spatial coherence be- 
tween the same regions both exist for an identical set of 
parameters. 



VI. DISCUSSION AND CONCLUSION 

Let us summarise the results of this paper. Initially we 
used a purity measure of entanglement to derive an en- 
tanglement witness, £, Eq. (|7]), that indicated whether 
two spatial modes of a non-interacting Bosonic gas of 
fixed particle number were entangled. The witness, £, 
is paramcterised by the one-body reduced density ma- 
trix, p{r,r'), that measures the spatial coherence be- 
tween points r and r ', and the probability to occupy 
one of the regions p. When the two regions occupy the 
entire confining geometry (p — I /2) the entanglement 
witness, £, depends only on the one-body density matrix 
and entanglement and spatial coherence are equivalent. 
If the probability to occupy a region is less than a half, 
p < 1/2, then spatial coherence between the regions is 
only necessary for entanglement. 

We plotted the results for the case of a cigar-shaped 
harmonic potential and as expected entanglement only 
exists below the critical temperature for condensation, 
Tc, when there is coherence between two separated re- 
gions. The amount of entanglement increases when the 
spatial modes are closer together because the population 
is located at the centre of the trap. If there were an equal 
probability to find Bosons at every point in the trap there 
would be an equal amount of entanglement between all 
regions of space. The non-local quantum part of the one- 
bodied reduced density matrix is as expected very long 
ranged. 

To gain a greater understanding of entanglement in the 
BEC phase transition it was also useful to consider multi- 
mode entanglement, as three or more spatial modes can 
still be entangled even when two spatial modes are not. 
We derived a second entanglement witness, W, to detect 
entanglement between three spatial modes using an ex- 
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isting entanglement witness, Tlw- In this case the three 
regions occupied the entire confining volume so the wit- 
ness, W was parameterised exclusively by the one-body 
reduced density matrix. The probability to occupy all 
three regions must be of the same order for entangle- 
ment to exist. At zero temperature, T = 0, the state 
with maximum entanglement corresponds to the region 
configuration that allows the probabilities for occupation 
of each region to be equal. Pa — Pb ^ Pc- The entan- 
glement in this state was also more robust to mixing and 
survived until the highest temperatures. 

In the introduction we spoke about two conjectures 
concerning the entanglement in quantum phase transi- 
tions. Namely that entanglement is responsible for the 
long-range correlations that develop in the transverse 
Ising and XY models at the quantum critical point and 
that at such critical points the system is maximally, or 
critically, entangled f2] . We will discuss these conjectures 
in the context of the BEC phase transition. Even though 
the phase transitions are very different in nature, our re- 
sults may act as a guide for how entanglement behaves 
in other such critical systems. 

In contrast to a quantum phase transition, the BEC 
phase transition is classical and is driven by thermal fiuc- 
tuations. One could, however, induce the BEC phase 
transition at zero temperature by increasing the density 
of Bosons until their wavefunctions start to overlap, and 
a BEC was formed. It is therefore sensible to compare the 
entanglement properties of BECs at zero temperature to 
entanglement in quantum phase transitions. Another dif- 
ference between entanglement in BECs and entanglement 
in discrete lattice type systems is the couplings between 
sites. In discrete systems one normally considers nearest 
neighbour couplings whereas in a BEC there are pairwise 
couplings between all regions. The entanglement struc- 
ture is far richer in a BEC, with entanglement stretching 
well beyond neighbouring regions. 

Entanglement below the critical temperature for con- 
densation occurs due to the long-range correlations that 
distribute Bosons coherently over the spatial modes. We 
have seen that spatial coherence is a precursor to en- 
tanglement and we know, of course, that the entangle- 
ment does not cause the spatial correlations. When the 



probability for Bosons to occupy each region is equal, 
entanglement and spatial coherence between the regions 
are equivalent and the system has the maximum possi- 
ble amount of entanglement. In discrete systems, where 
the sites are pre-defined, could entanglement and spatial 
coherence also be equivalent? At criticality long-range 
correlations do occur throughout the system, but one 
would need a suitable multi-partite entanglement wit- 
ness to fully characterise the behaviour of entanglement 
between all the sites at that point. 

Let us finally discuss two possible extensions of this 
work. Firstly one could characterise the behaviour of 
entanglement between M spatial modes. We know that 
genuine multi-mode entanglement only exists below the 
critical temperature for condensation, due to the require- 
ment of spatial coherence between distant regions. How- 
ever, it is interesting to ask whether there could be en- 
tanglement above the critical temperature for conden- 
sation. Consider regions smaller than the Boson's de 
Broglie wavelength, so that there was spatial coherence 
between neighbouring modes, entanglement of the whole 
system could then be built up from bipartite entangle- 
ment, to form some sort of entangled graph state. 

A second reason for understanding how spatial coher- 
ence is related to entanglement lies in the fact that it 
is becoming increasingly accepted that BEC like coher- 
ences can occur in systems at higher temperatures. In 
this case the systems interacts strongly with the environ- 
ment and a constant supply of energy lift the systems 
modes out of the Planck distribution causing the strong 
excitation of a single mode [TU] . This is particularly inter- 
esting when one realises that biological systems may show 
such behaviour where quantum coherence is thought to 
be present in energy transfer in photosynthetic systems 
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